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Abstract 



The fermionic, bosonic and supersymmetric variants of the colour-flavour transformation 
are derived for the orthogonal group. These transformations are then used to calculate the 
ensemble averages of characteristic polynomials of real random matrices. 

1 Statement of main results 

Since the pioneering work of Wigner, many physical systems have been successfully studied 
with the help of random matrix models. Among these asymmetric real random matrices 
arising in applications in neural networks [I], quantum chaos [2] and QCD [3] are known 
to be the most difficult. In fact, until the very recent breakthrough H @ U, the 
eigenvalue correlation functions of real and complex eigenvalues were not accessible even 
for Gaussian matrices and calculating the ensemble averages of eigenvalue statistics in the 
complex plane for a sufficiently general class of real random matrices remains a challenging 
problem. The mathematical difficulties in calculating the eigenvalue correlation functions 
of real random matrices are mainly due to the fact that their eigenvalues are either real 
or pairwise complex conjugate and the mathematical tools available to study real random 
matrices are very limited, especially when compared to those available for complex matrices. 
In this paper we derive several integral transformations dealing with integrations over real 
orthogonal matrices which we believe might be useful in the above context. These integral 
transformations are known under the name of the Colour-Flavour Transformations. 

The Colour-Flavor Transformations (CFT) are certain types of integral transformations 
based on Howe's dual pair theory. They were was first derived by Zirnbauer [9] in 1996 and 
since then have became a standard tool in mesoscopic physics, random matrix theory, lattice 
QCD as well as other fields. The CFTs were originally derived for V(N) and Sp(2iV), and 
later generalized to other classical groups [THl EH [HI [13] . However, both the fermionic and 
bosonic variants of the O(N) CFT appearing in the literature [HI [13] do not seem to be 
reproduced in the correct form, possibly suffering from typos. In this paper, we first correct 



the bosonic and fermionic versions of the O(N) CFT and then derive the supersymmctric 
version, which is a new result. 



Fermionic 0(N) CFT 

dO exp$*O y ^) = Cjf f d»(Z, tf) exp \ fez o6 $ + ^Z^l) , (1.1) 
O(JV) Jz=-Z T z \ J 

where Cq is the normalization constant dependent on N and n, see Eg. (|2.20p and 

dfx(Z) Z t) = dZdZ] . (L2) 

Integral on the left-hand side of Eq. (|l.ip is over the real orthogonal group O(N). Here, 
ip l a and ^> % a , i = 1, . . . , N and a — 1, . . . , n are Grassmann variables. Note that the V>'s are 
not necessarily related to ^'s by complex conjugation. That is, one can replace tp l a with an 
arbitrary set of independent Grassmann variables which are not related to "0 by any operation 
and the identity will still hold. The integral on right-hand side of Eq. (|l.l[) is over complex 
skew-symmetric matrices of dimension n x n. 

Bosonic 0(N) CFT 

dO eM&Oijti) = Of / MZ, ^) \ UaZa^l + fiZ^i) , (1.3) 

O(JV) Jl-ZZi>0 1 \ J 

where 

dfj,(Z, Z ] ) = det^^-^l - ZZ*)dZdZ* . (1.4) 

Integral on the left-hand side of Eq. (|1.3p is again over the real orthogonal group O(N). Here, 
<f) l a and <j) l a , i = 1, . . . , N and a — 1, . . . , n are complex variables. The normalization constant 
Cq is defined in Eq. (|2.26[ ). Similar to the fermionic case, one can also replace <f> l a with an 
arbitrary set of independent complex variables and the identity will still hold. The integral 
on the right-hand side of Eq. (jl.3p is over complex symmetric matrices of dimension n x n 
such that 1 — ZZ' is positive definite. 

Supersymmetric 0(N) CFT 

dOexp(COy^)= f d^Zj)cxpU^lZ ab ^l + ^l(-lpZ ab rb) ■ (1.5) 

O(N) JM B x.M F 1 \ / 

As before, the integration on the left-hand side is over O(N) and i,j — 1,...,N. Here 
ip and ip are graded vectors whose elements ip^ and ij) l a are bosonic when a = (a, B) and 
fermionic when a = (a, F), where a — 1, ... ,n. On the right-hand side, Z and Z are In x 2n 
dimensional supermatrices subject to the following condition 

Z = Z T a, Z = (tZ t , (1.6) 

where a is the superparity defined as a = ( ^ ^ V Here, Z and Z parameterize the 
Riemannian symmetric superspace CI|DIII defined as OSp(27i|2n)/GL(n|n), such that 

M B = Sp(2n,,M)/U(n), M F = SO(2n)/U(n) . (1.7) 
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We defined dfi(Z,Z) = dZdZ SDet T (l — ZZ), where dZdZ is the flat Berezin measure on 
the space of supermatrices Z and Z. As usual, SDet denotes the superdeterminant. The 
integration domain Mg x Mp is fixed by 

Zbb = Zg B , Zpp = —Z FF . (1-8) 

The symbol \a\ denotes the parity of index a, \a\ = when a is bosonic and \a\ = 1 when a 
is fermionic. 

With the fermionic CFT, we calculated the ensemble average of spectral determinants of 
certain real random matrices. The first one is the positive integer moments of modulus of the 
characteristic polynomial of matrices of the form GO (TV), where G is diagonal. Its average 
over the orthogonal group 0(N) can be written as an integral of products of Pfafnans, 



(|* - GO\ 2m ) 0{N) = [ dO det m [(z - GO){z - GO) 1 "] 
Join) 



'O(N) 

N I" r, — | 

(1.9) 



const. / dn(Z,Z*)f\pf 



Im. Z^ 



Here we introduced 2x2 matrix £t? — diag(z, z) and the integration measure is defined 
in Eq. (jl.2p . In the second example, we calculated the ensemble average of characteristic 
polynomials of real random matrices A from the Jacobi ensemble, 

I w {\i) = J dAdA T det a (AA T )det b (1 — AA T ) det(A — A) det(7 — A T ) , (1.10) 

where a and b are non- negative integers. This average is again written in terms of Pfaffians, 

f 1 1 

Iw(\l) = const. / dr pf [ai,j]i,j=o,...,2s-i , when N — 2s . (1-H) 

Jo I 1 + r ) 



and 



I W (\, 7 ) = const. ^ dr (i+r)jv+2 pf _ %(a ^. 1} Q 



, when N = 2s + 1 . 

ij'=0,...,2s 

(1.12) 



Here is defined in Eq. (|3.43[) and h is defined in Eq. (|3.42p . 

This paper is organized as follows. In section 2 we give a detailed proof of the fermionic 
CFT and then outline the derivations of the bosonic and supersymmetric CFT. In section 
3, we calculate the ensemble average of characteristic polynomials as applications of the 
fermionic CFT. A summary is given in section 4. 



2 Proof of colour-flavor transformation 

We derive the O(N) colour-flavor transformations in this section. A detailed derivation is 
given for the fermionic case. Derivations for the bosonic and supersymmetric cases are given 
in more brief way with only important issues addressed, whereas details are referred to either 
the fermionic case or literatures. 

In this paper, all the three types of colour-flavor transformations are established by 
Howe's dual pair theory, see [9l [15] and references therein. It is worth mentioning that in 
certain cases symmetric polynomials can also be used to derive these transformations or even 
transformations with different forms |11U16| . However, it turns out that here the derivations 
based on dual pair theory are more convenient. 



3 



2.1 Fermionic 0(N) CFT 

To establish the fermionic colour- flavor transformation over O(iV), we first derive the trans- 
formation over the special orthogonal group SO (TV), 



SO(N) 



d/i(0) exp(^Oy^) = C [ dfi(Z, Zt) exp ± fez o6 $ + (1 + ^det X) , 

(2.13) 



F C 

where CK and K = „ F * ... are constants defined later and M. is an N x N matrix defined 

°0 A l iV! 



as Mij = 



r a (i + zz% b ^ b 



Remark: To get Eq. dTTTJ) , we exploit the fact that 0(N) = + (N)®0-(N), where 0+ = 
SO (AT) and 0_ is a rotation followed by a 'reflection' i?, which can be chosen as R = 
diag(Ijv-i) —1). Note that 4i l Rij flops the sign of all ip N, s therefore inverts the sign of det Ai. 
Hence, the parts containing det Ai are cancelled when we combine the contributions from 
the normal rotation SO (AT) and the improper rotation R ■ SO(Af). Which proofs Eq. fll.ljl . 
In fact, it turns out that one can 'naively' use a similar mapping as defined above Eq. (|2.14p 
for 0(N) and get Eq. (|l.l|) with less effort. However, due to the fact that detO = ±1 for 
O € O(AT) this kind of mapping is mathematically better defined for SO(N). 

In the following paragraphs, we will mainly follow the method of paper [lOj . Introduce 
fermionic creation and annihilation operator /* and /*, where i = 1, . . . , N and a = 1, . . . , n. 
As usual, we borrow the terminology from lattice gauge theory where the upper indices 
are referred to as 'colour' and the lower indices are referred to as 'flavor'. This set of 
operators satisfy the canonical fermion anticommutation relations = & lJ 5ab and 

{/*, fl} = {fa, fl} = 0. These operators construct a Fock space &p for a fermionic system. 
Let |0) be the vacuum state, then we have /q|0) = 0. The quadratic operators f a f\, flfl and 
fafl ~ fafh therefore define a representation of the Lie algebra so(2nN, C). This algebra has 
to two commuting subalgebras. Which are so(2n, C) generated by flfl, flfl and p a fl — f a fl 
and so{N, C) generated by ft ft + Pjl- 

To derive colour-flavor transformation we need to construct two types of projection op- 
erators Pc and Pp which project states to a subspace of &p, named colour-single space, 
which is invariant under the SO (AT) group action O i— » To '■= exp(/*(ln 0)ijfl). As usual, 
we define 

dO T ■ (2.14) 

SO(AT) 

Note that for SO (AT), the colour-single space has two disconnected components, each of 
which carries an irreducible representation of SO(AT). One of them contains the vacuum 
state |-0o) = |0) and is in addition O(N) invariant, and the other one contains the state 
\ipi) :— flff . . ./i^O)- Acting on the vacuum and \ipi) by the operators span each 
colour-single subspace. Therefore, we have Pp = Pp + Pp 1 . 
It is convenient to define the notation 

^ _ / f\ A = l,...,n , j _ j f_\ A = l,...,n 

A ~{r A - n A = n+l,...,2n &nd ^ ~ \ f A _ n A = n+l,...,2n ' (2 ' 15) 

For geSO(2n), It is direct to check the mapping g i— > T g := exp(^c l A (lng)ABC l B ) constructs 
a representation of SO(2n), such that T g S , A T' g 1 — c^gsA- An element of SO(2n) can be 
parameterized as [18] 

9 = ( ? where uu ^ + vy] = 7 ' and uyT + vljT = ' ( - 2 ' 16 ' ) 
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Note that this G is a subgroup of U(2n) that is isomorphic to SO(2n). It is clear that G has 
a U(n) subgroup G D H = diag(t/, U) = U(n). We parametrize the coset space G/H as 

G,H ~ + ztf)-l (i + tfz)-i ) ' (2 - 17) 

where Z = U~ 1 V are complex skew-symmetric matrices of dimension n x n. 

p F = E p r° = E C ° I d 9 TglfaK^T- 1 , (2.18) 

<x=0,l o=0,l - , SO(2n) 

where is defined by the normalization condition (ipcr\PF a l^/v) = !• Define Z) = 
exp(i/*Z a h/^)|0) then it is direct to check 

Pf = C F / d ff T^lVo)^^" 1 = Co / d/i(Z, tf) |Z)(Z| . (2.19) 

JSO(2n) JZ=-Z T 

Using the method in chapter 2 of [17], we get the normalization factor 

c "=- 11 r(iv + i) • (2 ' 20) 

Note that we choose the group volume of SO(iV) to be 1. Follow the method in 10], C-f 1 
equals to the dimension of the representation of SO(2n) defined by the Young diagram with 
1 row and N columns. 

Now define — exp(/*?/;* )|0), then the colour-flavor transformation follows from the 
identity (ip\Pc\ij>) = (iP\Pf\iP)- It is straightforward to show 



&\P c \i>)= dOeM^lOi^i), (2.21) 

^SO(N) 



and 



(^\P Fo \xl>) = Cjf J ^ d^Z, tf) exp 1 (fcZabft + tizlti) ■ (2-22) 

Follow the method in [10] again, we get 

(4>\P Fl W) = ^ f d^{Z, Zt) exp i fez o6 $ + ^Z^) det [#(1 + ZZ^U^ 

(2.23) 

where A" is the dimension of the representation of U(n) defined by the Young diagram with 
1 row and N columns. Combining the above formula, we complete the proof of Eq. (|2.13[) 

2.2 Bosonic O(N) CFT 

The bosonic colour-flavor transformation is derived with almost the same method |13] but 
with the following few changes. We construct a bosonic Fock space &b with bosonic opera- 
tors. The two commuting algebras are so(iV, C) and sp(2n, C) [13]. Define 

b\ A = l,...,n , j _ / b\ A = l,...,n 



c l A = { ,i a o and C A = { \i a Ti' o ■ (2-24) 

\ b\_ n A = n+l,...,2n A \ -b\_ n A = n + l,...,2n y ' 
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For <?GSp(2n, R), it is straightforward to check that the mapping g i— > T g := exp(ic^(ln g)ABC B ) 
constructs a representation of Sp(2n,R), such that T^c^T" 1 = c B gsA- An element of 
Sp(2n, R) can be parameterized as 



U V 
V u 



where UU f - V T V = I, and U f V = V T U . (2.25) 



Here G is a subgroup of U(n,n) that is isomorphic to Sp(2n,R). Note that G has a U(n) 
subgroup. The coherent states are similarly defined as before, \Z) = cxp(^b l a Z a i,bl)\0) and 
\4>) = exp(6* <j>D 1 0). Following the same procedure as before and [13], we get Eq. (|1.3p . 

In order to fix the normalization factor, choosing the group volume of O(N) to be 1 and 
using Eq. of [T7], we get 

_^ 1 ) iV-2n " rr 1 (f -j)T(N-l-i) 

c " =n —2~ 11 r(jv-i-2i) ' (2 ' 26) 
2.3 Supersymmetric O(N) CFT 

When -0 and ip are graded vectors, we need to extend the flavor group to incorporate this 
'supersymmctry'. The Lie superalgebra we need satisfies the following condition 

Q = -lQ T l ■ (2-27) 
Let cr's being the Pauli matrices, 7 can be defined as 

7 =(^ °)®I W . (2.28) 

Therefore, Q defines the Lie superalgebra osp(2n|2n) whose boson-boson block is the Lie 
algebra sp(2n) and fermi-fermi block is so(2n). The remaining derivation closely follows [9]. 
To do a simple check, setting the bosonic (or fermionic) degree of freedom on both sides of 
Eq. (|1.5p to zero and choose corresponding integral measure, we can recover the fermionic (or 
bosonic) CFT Eq.{THJ. 



3 Applications in real random matrices 

In this section, we calculate the averaged characteristic polynomials of real random matrices 
from certain ensembles. These calculations are simple applications of the fermionic O(N) 
CFT. 

3.1 Modulus square of characteristic polynomials of certain type of 
real random matrices averaged over O(N) 

Let z be a complex number and G be a real diagonal matrix. We calculate the following 
quantity, 

F G (z) = (\z-GO\ 2m ) 0[N) = f dO det m [(2- GO){z-GOy\ . (3.29) 

JO(N) 

Note that (\z — GU\ 2m ) V ( N ^ where the average is over unitary group, has been calculated 
with the method of symmetric functions in [21] . And it is not hard to show that the method 
we will use in this section can be applied to the unitary case as well, in which case one will 
use the XJ(N) fermionic CFT. 
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Introducing Grassmann vectors, tf a and i — 1, . . . , N and a — I, . . . ,m, we can re- write 
Eq. (ET34l) as 



F G (z) = / dO [ dfjdr]d^ e - z ^-^^ Bi0i ^ + ^ *i3iii 

Jo(N) J 

= [ dfjdr,dtdti e"*«-^ / dDexp^jO^), (3.30) 

J JO(N) 

where we introduced the composite notation {4>\ a: V>| a ) = (Va9ii ~C a 9i) ( n0 summation) and 
{ip\ a ,tp2a) = (Va'Ca) an d changed the order of integration. Defining the 2m x 2m complex 
skew-symmetric matrix Z and applying the fermionic CFT Eq. (|l.ip , we get 

F G (z) = Qf J dfjdrjd^ e-**"'-*^ ^ ^ t d/i(Z, Zt) exp ± fez o6 $ + ^Z^j) ■ 

(3.31) 

Next, we change the order of integration and use the standard Gaussian integral formula for 
Grassmann vectors. The final result is written in terms of Pfaffians and given in Eq. (|1.9|) . 
When ttl — 1 , we can integrate over Z explicitly since 



a 

v -a 

By Eq. (|3"3Tj) . we get 

1 N 

F G {z) = (\z -GO\ 2 )o(N) = const. / dada + a -, N+2 Y[( 



(3.32) 



z\ 2 + aag 2 ) 



cons t. ■£r\z\* lN - l) S l (Gf 2 ) , (3.33) 

Here, C l N = N\/l\(N — 1)1 is the binomial function and S l (G 2 ) denotes the 1-th order elemen- 
tary symmetric polynomials of g$% i.e. S°(G 2 ) = 1, S^G 2 ) = S 2 (G 2 ) = T,i<j 9^1 



etc. 



3.2 Characteristic polynomials averaged over Jacobi ensemble 

As another simple application of the transformations we derived in the previous section, we 
calculate the following average, 



/ 



dAdA T W(AA T ) det(A - A) det( 7 - A T ) . (3.34) 



Where W(AA T ) is an arbitrary invariant ensemble which depends on singular values of A 
only, i.e. W(A) = W(OAO') = W(G), where O and O' are arbitrary elements in 0(N) and 
G are the singular values of A defined by A = 0\G02- In particular, we are interested in 
ensembles whose potentials are also separable functions of G, i.e. 

W(AA T ) = Y\_W{g 2 ) ■ (3.35) 

i 

First, we make singular decomposition A — O1GO2, where G = diag(<?i, . . . , (/at) > 0. 

j dAdA T W{AA T ) = j dO x d0 2 J H W(g 2 )d 9i [] \ g * - g 2 \ . (3.36) 



i<3 
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By invariance of Haar measure, we have det(A — O1GO2) = det(A — O2O1G) = det(A — OG). 
Introducing Grassmann vectors as before, we can re- write Eq. (|3.34j) as 



O(N) 



dO J J J II </</, J [ g 1 g] j df)dr]d^di, e A " " " " " " 1 " 9 :; . 

(3.37) 



Use the same method as in Eq. (|3.30p -Eq. (|3.31[) . we get for separable invariant potentials 
Eq.dSSSJ), 



J W (A,7)= J dAdA T W(AA T ) det(A — A) det(7 — A T ) 

=const. I' dr ^ / ft Iff? - 3j\ U (*7 + ^X) W(gf)dg t . (3.38) 

' i<j i 

One benefit of using singular values instead of eigenvalues is that all integral variables are 
real non-negative. From the above formula it is clear that when 7 = A the integral on the 
left-hand side depends only on |A| 2 . 

Note that Eq. (|3.38[) applies to all ensembles of real random matrices satisfying Eq. (|3.35p . 
These ensembles include W(AA T ) — e~ TlV ( AA ' T \ which becomes the Ginibre ensemble when 
V(x) = ix, which has been studied intensively. In the remaining part of this section, we 
consider the Jacobi ensemble, 

W{x) =x a {l-x) b . (3.39) 
To proceed, we use the method in Chapter 5 of [19] . 

f \[g: //J Y[(l + rgf)W(g 2 )d gi = Nl J Y[d 9i det [^( 5 ?)^_x( 5l 2 )(l + r^)] . 

i<3 i 0< SN ---<9l<l i 

(3.40) 

where Rj (x) can be arbitrary monic j order polynomials of x. Here we choose Rj (x) = x J . 
Define the following functions 

h(a,b;x)= dgg^{l-g^ = -T{b+l)T{a+-)Y J ^ : 7 : r , J^a, > ^Al) 
Jo 2 2 j-^ T(b - 1 + l)T(a + 1 + %) 



and 



h(a,b;x)= / dg {1 + rg 2 )g 2(a+l) (1 - g 2 ) b = h{a + i,b; x) + — h(a + i + l,b;x) , (3.42) 
Jo A 7 

Let B(x, y) — be the Beta function. When N = 2s, introduce the skew-symmetric 
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matrix a, 



<*ij = / dg (1 + rg 2 )g 2a (l - g 2 ) h {g*%{g) - g^h(g)) 
Jo 

_1 f B(6 + l,a + j + |) 5(26 -Z + l,2a + l+ i+ j + i) 
_ 4\ a + 6 + j + | ^ B^-Z + M + j + Z + f) 

r B(b+l,a + j + |) A B(26 - ? + 1, 2a + 2 + i + j + Q 
A7 a + 6 + j + § ^ B(6-/ + l,a + j + / + §) 

r 5(6 + 1, a + j + |) A B(26 - i + 1, 2a + 2 + i + j + Q 
A7 a + b + j + l B(6-/ + l,a + j + / + §) 

/ r A 2 B(b+ l,a + j + |) A 5(26- 1 + l,2o + 3 + i + j + I) . .1 
Ia^J a + 6 + j + f ^ B(6-Z + l,a + j + / + §) z ^/ ' 

(3.43) 



By Eq. (5.5.8) and Eq. (5.5.9) of [32], we get Eq.JTn} and Eq.JTEJl. 

Similarly, for W{x) = exp(— we can check that Eq. (|3.38p and Eq. (|3.40p gives the 
same results in [30] 

f dAdA T e-? TlAAT dct(A - A) dct( 7 - A T ) = const. ^ . (3.44) 



4 Conclusions 

In this paper, we derived the 0(N) colour-flavor transformations. We believe that these 
transformations will be useful in the study of real random matrices, just as the U(iV) CFT are 
for complex random matrices. As a simple application, we calculated averaged characteristic 
polynomials for two types of ensembles, which, to the best of our knowledge, have not been 
calculated before. 

We have only showed examples with the fermionic CFT whereas applications of the 
bosonic and supersymmetric CFTs will be postponed to future works. In which, as usual, 
we need to pay attention to analytical issues, especially when we need to change the order 
of integrations, [211 [22]. Also, by its nature the supersymmetric CFT often requires more 
work. However, when n = 1 Eq. (|1.5[) enjoys the simplicity that 

*-U °) "° d <>)■ <" 5) 

This is because 0(2)/U(l) consists of only single points. 
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